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ABSTRACT 

Archival lUE spectra of the VY Scl system MV Lyr, taken during an inter- 
mediate state, can be best fit by an isothermal accretion disk extending half way 
to the tidal cutoff radius. In contrast, a recent HST spectrum, while MV Lyr 
was in a high state, can be best fit with a standard T{R) profile for an accretion 
disk extending from an inner truncation radius to an intermediate radius with 
an isothermal accretion disk beyond. These fits use component star parameters 
determined from a study of MV Lyr in a low state. Model systems containing 



^ Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc. under NASA contract NAS 5-26555 and with the Apache Point Observatory 3.5 m telescope, which is 
operated by the Astrophysical Research Corporation. 
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accretion disks with standard T{R) profiles have continua that are too blue. The 
observed high state absorption line spectrum exhibits excitation higher than pro- 
vided by the T{R) profile, indicating likely line formation in a high temperature 
region extending vertically above the accretion disk. The absorption lines show a 
blue shift and line broadening corresponding to formation in a low velocity wind 
apparently coextensive with the high temperature region. Lines of N V, Si IV, C 
IV, and He II are anomalously strong relative to our synthetic spectra, indicating 
possible composition effects, but unmodeled excitation effects could also produce 
the anomalies. 

An analysis of a low state of MV Lyrae, considered in an earlier study and 
extended in this paper, sets a limit of 2500K for the Teg of an accretion disk that 
may be present in the low state. This limit is in conflict with two recent models 
of the VY Scl phenomenon. 

Subject headings: accretion, accretion disks — novae, cataclysmic variables — 
stars: individual(MV Lyrae) — ultraviolet: stars — white dwarfs 



Introduction 



Cataclysmic variables (CVs) are binary stars consisting of a white dwarf (WD) and a 
late spectral type secondary which fills, or almost fil ls, its Roche lo be and transfers mass 
to the WD. CVs divide into a number of subclasses (j Warner 1 1 19951 ) and the object of our 
study, MV Lyrae, is a novalike CV of the VY Scl subclass. A defining characteristic of this 
subclass is a non-periodic alternation between high and low luminosity states, differentiated 
by a change in mean brightness of more than one magnitude between states. This alternation 
usually takes place on intervals of hundreds of days. 

The orbital pe riod of MV Lyr, 3.19'^, places it close to the long-per iod end of the period 



gap ( 


Warner 


1995) 


(2004 


). While it is 



gap is controlled by magnetic braking (jPatterson. 



of the gap is driven by gravitational radiation (IKing 



198 



King Ill988l) , and evolution shortward 



19881) . this scenario has b e en re cently 



challenged by, e.g ., Andronov. Pinsonneault. fc Sills (j2003 ). King fc Schenker (2002), and 
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Schenker fc King I ( 20021 ). Patterson (1984) developed an empirical mass transfer relationship 
for CVs which predicts a mass transfer rate of 2.5x 10~"'^°A^(^ yr~^ at the orbital period of 
MV Lyr. A theory of magnetic braking (IHameury et al. Ill988l ) shows that the mass transfer 
rate from the secondary is a decreasing function with decreasing orbital period, reaching a 
value of 1.0xlO~^A^oyr^^ at the top of the period gap, where MV Lyr is found, for the 
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parameters used in their study. Some authors, as discussed below, adopt even larger mass 
transfer rates at the top of the period gap. 

If the magnetic braking mechanism that drives evolution of the secondary star is sud- 
denly turned off at an orbital period of 3*^, then the secondary approaches thermal equilibrium 
on a Kelvin -Helmholtz ti mescale. Consequently, mass transfer continues at a decreasing rate 
for ^lO^yr ( iRitter Ill988l ). But the alternation between high and low states in VY Scl stars is 
on a much s horter timescale a nd re quires a separate mechanism to switch the mass transfer 
on and off. iLivio &: Pringle I (1l994j ) discuss the problem and propose starspots on the sec- 
ondary star, which pass acro ss the LI point and caus e reduced mass transfer. This proposal 
has been further developed (IKing &: Cannizzo 1119981) with application of a time-dependen t 
code to model the variable a ccreti on disk. (See also lSchreiber. Gansicke. fc Cannizzo Il2000l ). 
But the King fc Cannizzo I (119981 ) model produced the unwanted result of dwarf-nova type 
outbursts in the low state, after mass transfer had been switched off. A suggested escape 
was maintenance of the inner accretion disk in a permanent high state through irradiation 
by the WD (Leach et al. 1999, hereafter L1999). L1999 concluded that the presence of a 
40,000K WD does suppress dwarf nova outbursts. The L1999 study adopted a high state 
mass transfer rate of l.lxlO~*Al0yr~^ that was chosen to guarantee that hydrogen is fully 
ionized at the outer accretion disk boundary. According to L1999, cessation of mass transfer 
produces a rapid transition to a state of low but non-zero viscosity. The accretion disk re- 
mains nearly intact during the low state, transferring only a few percent of the accretion disk 
mass to the WD. Because of non-zero viscosity, the accretion disk still radiates at a low level 
during the low state. Hameury & Lasota (2002, hereafter HL02) show that the adopted WD 
mass in L1999 (O.4A^0) is important to the suppression of outbursts, and that a WD mass 
of O.7A^0 would have produced outbursts. HL02 alternatively propose that VY Scl stars 
have magnetic WDs and that the magnetic field truncates the accretion disk. HL02 show 
that, for a WD mass of 0.7 A4q and a magnetic moment of /i = SxlO^'^Gcm^, no outbursts 
occur and the change in the V magnitude, from high state to low state, is 5 magnitudes. 



Hoard et al. (2004, hereafter H2004) used the BINSYN suite flLinnell fc Hubenv Ill996f ) 
to calculate a system model and corresponding synthetic spectrum for MV Lyr in a recent low 
state. The synthetic spectrum accurately fits FUSE spectra and contemporaneous optical 
spectra, as well as lUE spectra from a prior low state. The H2004 discussion considers the 
connection of the low state study to earlier investigations of MV Lyr and shows that the 
low state can be understood in terms of a naked hot WD with a temperature of 47,000K, 
a photosphere log g of 8.25, and a metallicity of 0.3 solar. T he log q, together with the 
mass-radius relation for zero-temperature carbon WD models (IHamada fc Salpeter Ill96ll ) 
selects a model with A^wd = O.73A^0, estimated error of O.1A^0, and radius of O.OlO67i?0. 
The secondary star fills its Roche lobe and is cooler than 3500K; it contributes nothing to 
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the FUV flux and a varying amount to the optical flux, from 10% at 5200A to 60% at 7800A. 
If an accretion disk is present it contributes negligibly to the observed spectra. Assuming 
no disk is present, the model optical light curve, for a system orbital inclination of 12° taken 
from the literature, has an amplitude approximately 50% larger that of the observed optical 
light curve. The scaling of the model spectrum to the observed data leads to a distance of 
d = 505±50pc to MV Lyr. 

In §2 of this paper we describe the program suite used to analyze the observed spectra. 
In §3 we further investigate the low state and derive improved system parameters. In §4 we 
show that it is possible to augment the low state system with an accretion disk with artificially 
elevated T^s at large radii and achieve a close fit to intermediate state lUE archival spectra. 
In §5 we show that our models produce a predicted AV, low state to intermediate state, 
that is one magnitude larger than observed, indicating a problem with the models. In §6 we 
show that an isothermal 14,000K accretion disk extending half way to the tidal cutoff radius 
resolves the AV problem and provides a fair fit to the lUE spectra. §7 presents a model fit 
to a new high state HST spectrum, §8 considers the predicted luminosity change low state 
to high state, a general discussion is in §9 and §10 summarizes our results. 



2. Description of the BINSYN program suite 



The BINSYN suite (ILinnell fc Hubeny Ill996l ). used in the analysis presented in this 



paper, is a set of programs for calculating synthetic spectra and synthetic light curves of 
binary star systems with or without accretion disks and in circular or eccentric orbits. The 
model stars have photospheres of assigned Roche potentials, defined by grids of points with 
associated photospheric segments. Program modules assign T^s values to individual stellar 
points based on adopted polar Tes values, bolometric albedos, and gravity darkening expo- 
nents. Model stars may rotate synchronously or non-synchronously up to critical rotation. 
Program modules calculate log g values at stellar photosphere points based on adopted Roche 
potentials and stellar masses, while other modules determine whether a given point is visible 
to the observer, for a point on a star, the accretion disk face or the accretion disk rim, and 
for the current values of orbital inclination and orbital longitude. Calculation of a synthetic 
spectrum uses a set of precalculated source synthetic spectra for each system object (stars, 
accretion disk face, accretion disk rim). The source spectra for a given object typically form 
a two-dimensional set in Teg and log g that encompasses the extreme T^s and log g values on 
the photosphere of the given object. Program modules interpolate an individual synthetic 
spectrum to each stellar photospheric point and integrate over that star, with allowance for 
point visibility; the routines properly allow for Doppler shift and line broadening. 
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An accretion disk is specified in BINSYN by an inner and outer radius, and the accretion 
disk is divided into a specified number of cylindrical annuli, with each annulus divided into 
a specified number of segments. Radiation characteristics of the accretion disk depend on 
a temperature profile, T{R), of the accretion disk annuli. We reserve the term "Standard 
Model" to designate a T{R) relation 



T{R) 



with /3 = —0.75 for the Standard Model (IFrank. King, fc Raine Ill995l . hereafter FKR). W is 
a normalizing factor, (not a dilution factor); W = 1.0 for the Standard Model. For default 
conditions, our definition is algebraically identical to the FKR definition. 

The accretion disk is flared, in the sense that successive annuli step up in height between 
the inner radius and a specified semi-height at the outer radius. Let H be the semi-thickness 
at the outer radius r^, and let Vh be the inner radius. Then the semi-thickness at radius r, 
h{r), is h{r)/H = (r — rb)/{ra — r^). See Table 1. In calculating a model for an accretion 
disk system, BINSYN, by default, assigns Tes values to annuli based on the Standard Model, 
appropriate to an assumed mass transfer rate. An option permits modification of the assigned 
annulus T^fi values to include calculated irradiation by the WD, based on a bolometric 
albedo formalism. Calculation of a synthetic spectrum for the accretion disk, as part of 
the calculation of the system synthetic spectrum, involves interpolation among previously 
calculated source spectra (with their individual T^s values) to produce a synthetic spectrum 
for each annulus, for its T^g as assigned by BINSYN. The accretion disk sourc e spectra 



may c onsist largely of synthe tic spectra for annulus models produced by TLUSTY (IHubeny 



19881 : iHubeny fc Lanz I Il995l ). This calculation is followed by integration over all annuli, 
with detailed allowance for Doppler shifts and sources of line broadening, and with proper 
allowance for visibility to the observer. The number of source spectra typically will differ 
from the specified number of annuli. The programs finally output the calculated fiux as 
function of wavelength for the individual objects as well as their sum, which represents the 
system synthetic spectrum at the current orbital longitude. A parallel procedure produces 



synthetic light values, repeated for a grid of orbit al longitudes and so 



for specified wavelengths. Additional details are in iLinnell &: Hubeny I (119961 ) 



eads to light curves 



3. Fit to the low state optical spectra 

In fitting the optical spectra, H2004 found two difficulties: (1) The polar Teg of the 
secondary was cooler than 3500K, which is the temperature of the coolest (Kurucz) model 
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1500 2000 2500 3000 
Wavelength (A) 

Fig. 1. — Fit of synthetic WD spectrum to low state /[/i? spectrum of MV Lyr. The hght hue 
is the lUE spectrum, the heavy Une (colored blue in the electronic edition) is the synthetic 
spectrum. (See the electronic edition of the Journal for a color version of this figure.) 
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atmosphere available to us. (2) The optical spectra were observed through thin clouds so an 
empirical scaling factor was applied to the observed fluxes to match them to the synthetic 
spectrum (the latter was scaled to fit the FUV spectrum). 

Fig. 8 of H2004 shows the model fit to observed spectra, covering the interval 930A to 
7800A, including archived /C/E spectra SWP10905, covering 1150A to 1978A, and LWR09590, 
covering 1851A to 3349A, which were obtained on 1980 Dec. 27. Fig. 1 of this paper presents 
the model fit to the same low state lUE spectra in much greater detail than in H2004. Note 
that Fig. 8 of H2004 is a logaritmic plot while Fig. 1 of this paper plots fiux directly. In 
producing Fig. 1, the spectra were divided by the same scaling factors applied in the H2004 
plots (2.43x10^^ for the WD synthetic spectrum, 1.0x10"^^ for the /t^i? spectra), set by the 
fit to the FUSE spectrum. The mean residual of the fit to the 877 points between 1150A and 
3000A was 0.012 and the mean absolute residual was 0.038. These residuals apply to the 
scaled spectra. An important conclusion from the Fig. 1 fit is that there is no evidence for 
the presence of an accretion disk. In a test of the sensitivity of the fit to added fiux from 
a putative accretion disk, we find empirically that addition of 0.005 fiux units (Fig. 1 ordi- 
nates) produces a visually detectable displacement of the summed spectrum from the Fig. 1 
fit. At 3000A, the secondary star makes a calculated contribution of 0.1% to the system 
fiux. 

H2004 fit an optical (unfiltered) light curve of MV Lyr (their Fig. 7). This fit was derived 
from the BINSYN spectral model described above. It as s umes an orbital inclination for the 



CV ofi = 12° from lSkillman. Patterson, fc Thorstensen I (119951 . hereafter SPT95). However 



the corresponding amplitude of the model light curve is 50% too large. The orbital period 
light variation results from the variable presentation of the irradiated secondary component 
to the observer. As the assumed inclinat ion decreases from i = 12°, the light amplitude 



decreases monotonically. An improved fit ( ILinnell et al. 1120051 ) implies a system inclination 
of i = 7±1°, which we adopt for this analysis. 

A remark is in order on the absolute accuracy of the results obtained in this paper. A 
key step is in the determination of the WD mass and radius in units of the solar values by 
comparison between the MV Lyr WD log q f r om o ur analysis of the FUSE spectrum and 



the log g values for the iHamada fc Salpeter I (jl96lh carbon WD models. The mass ratio 



q, comes from relatively uncertain radial velocity measurements (SPT95), and this value 
determines the mass of the secondary star and, since the secondary fills its Roche lobe, the 
size of the secondary star. The fit of our synthetic spectra to FUV data, in all cases, applies 
a scaling divisor, 2.43x10^^, to observed spectra which were divided by 1.0x10"^^ (see the 
ordinate labels). The scaling divisor may be inaccurate by 10% to 15%. 



We present here an improved representation of the secondary star, using the NextGen 
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synthetic spectra for M stars, version 5 (iHauschildt et al. Ill999l ). We find a good fit to the 



observed spectra with a polar Teg = 2600K. Irradiation by the 47,000K WD produces a 
calculated LI point Tgfj = 5344K on the secondary star. BINSYN interpolated among seven 
NextGen synthetic spectra with Tgfj values between 2600K and 5400K to determine flux 
values for individual secondary star photospheric segments. The fit of the system synthetic 
spectrum to the optical spectra is shown in Fig. 2. The optical spectra were multiplied by 
an empirically determined factor of 2 to correct for obscuration by clouds. 

Recapitulating, we have developed a model for the low state of MV Lyr consisting of 
a 47,000K WD and a secondary star with a polar T^g of 2600K. The secondary star fills 
its Roche lobe and is irradiated by the WD. The model has no accretion disk, and differs 
from the model in H2004 only in the improved representation of the secondary star in the 
calculated synthetic spectra. This model produces an excellent fit to the FUSE data, the 
low state lUE data, and the optical data. We proceed to compare our model to the low state 
predictions of L1999 and HL02. 

L1999 (their Fig. 3) predict an accretion disk luminosity in the low state of F5500 = 
7.0x10^"^ erg s^^A^^. Fig. 3 plots the residual of the L1999 accretion disk (asterisk), for the 
505pc (H2004) distance of MV Lyr. The residual is F^^oo/d^, where d = 505pc. The L1999 
accretion disk flux at 5500A would have to be approximately 60 times fainter to accord with 
the residual plot. 

We calculated a synthetic spectrum for the example in HL02, with fi-^Q = 5, with 
-'^wd = 0.73x7V1q and mass ratio q = 0.43. The calculated inner radius for the HL02 
example (equation 4) is nn = 5.7x10^ cm and the MV Lyr WD radius is r^d = 7.9x10* cm. 
Consequently we set the inner accretion disk radius equal to the WD radius, the outer radius 
equal to the tidal truncation radius, and the mass transfer rate equal to = 2xl0^^gm s~^ 
(HL02 value). We used the synthetic spectrum of the calculated accretion disk alone in the 
following comparison. We divided the output flux in the synthetic spectrum by 100.0, the 
factor by which the HL02 model luminosity drops between the high and low states, to provide 
the HL02 model flux to compare with our low state data. As with L1999, this calculated low 
state luminosity is a residual, [(system model)- (BINSYN low state model)], shown in Fig. 3 
by a diamond. It could be argued that division by 100.0 for all wavelengths is too simplistic 
in transforming the HL02 model to an equivalent low state, particularly since we have no 
observational data in the 5000A to 6OOOA interval. Our Fig. 2 shows that the low state 
synthetic spectrum varies smoothly over this wavelength interval, as does our intermediate 
state model and high state model, discussed in the following sections (not shown in separate 
flgures). Thus, we believe the division by 100.0 is appropriate. The HL02 flux needs to be 
smaller by approximately a factor 5-10 to agree with the spectral residuals. Thus, we flnd a 
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4000 5000 6000 7000 

Wavelength (A) 

Fig. 2. — Fit of synthetic spectrum to low state optical spectra of MV Lyr. The light line 
indicates the optical spectra, the heavy line (colored blue in electronic edition) is the system 
synthetic spectrum. The secondary star contribution to the system synthetic spectrum is 
the lower spectrum in the plot. (See the electronic edition of the Journal for a color version 
of this figure.) 
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Fig. 3. — The 0(optical spectra)-C (synthetic spectrum) residual spectrum from Fig. 2. The 
asterisk marks the 5500A flux of the L1999 model (see text) if it is placed at the distance of 
MV Lyr, while the diamond marks the corresponding datum for the HL02 model (see text). 
The feature at 7600A is due to the molecular oxygen " A" absorption band in the terrestrial 
atmosphere. The emission lines are produced on the secondary star by irradiation from the 
primary star. 
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discrepancy between the predictions of both L1999 and HL02 for the MV Lyr low state and 
our fit to available low state data. 

How cool would an MV Lyr accretion disk have to be, in the low state, to avoid detection 
in optical spectra? Tests with an isothermal accretion disk, represented by a black body and 
truncated at r = 1.7xr„d (see the following text for a discussion of the truncation radius), 
demonstrate that an accretion disk must have a T^g less than 2500K to avoid a conflict with 
our optical spectra (a 2500K accret ion disk would provid e a detectable 3% of the system 
flux at 7800A). A study of TT Ari ( Gansicke et al. Ill999l ) obtained comparable results. A 
possible accretion disk would need to be truncated at r = 12r„d and have a T^s less than 
3000K to satisfy the low state optical spectra. 



4. Fit to the intermediate state lUE spectra 

The lUE archives include 10 spectra of MV Lyr. None of these were obtained during a 
high state, and two, SWP07296 and LWR06288, were obtained on Dec. 2, 1979 during an 
intermediate state. Intermediate states are roughly one magnitude fainter than high states. 

Fig. 8 of H2004 shows the intermediate state lUE spectra on the same plot as the 
fit to the FUSE//t/i?/optical low state spectra, for comparison. A striking feature of the 
intermediate state lUE spectra is their much lower (i.e., less blue) spectral gradient than the 
low state lUE spectra. This is surely due to the presence of an optically thick accretion disk. 



We calculated a s eries of annulus models, using TLUSTY version 200 (IHubeny Ill988 



Hubeny fc Lanz Ill995l ) for individual mass transfer rates of l.OxlO~^A^0yr~^, 4.OxlO~^A^0yr" 
and l.Ox lO~^Al0yr~^. These latter models were spaced in values of r/r„d between 2.5 and 
an outer value of r/r„d = 28.5, Teff=10,620K. TLUSTY includes a large number of op- 
tional control flags that control the program performance. We preserved the default control 
flags, and thereby produced LTE thin disk models with thin photospheres. Convergence 
problems intervened for la rger values of r/r^j. We used a viscosity parameter a = 0.1 



( IShakura fc Sunyaev Ill973l ) for all models. The models included opacity contributions from 



the 30 most abundant periodic table elements. We use d these models to calculate sourc e 



(see above) synthetic spectra, via SYNSPEC version 48 (IHubeny. Stefl fc Harmanec Ill985l ). 
for all models, with the same opacities used for the model calculations. We included syn- 
thetic spectra for Kurucz stellar model atmospheres, calculated with SYNSPEC, to handle 
assigned temperatures cooler than the 10,897K limit described above. Comparison among 
the annulus source synthetic spectra for the different tabular mass transfer rates showed that 
spectra for different mass transfer rates but similar T^g values are essentially identical. Con- 
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sequently, it is permissible to fix the source synthetic spectra of the annuh at some specific 
tabular mass transfer rate, but change the BINSYN mass transfer rate for a new simulation, 
and still produce a valid system synthetic spectrum. 

The accuracy of the accretion disk synthetic spectrum is sensitive to the Tes spacing 
of the source synthetic spectra. Linear interpolation between source spectra differing in Tes 
values by several thousand Kelvins can produce an inaccurate interpolated spectrum. After 
extensive experiment, we adopted a "universal" set of 69 source spectra, of which 62 were 
TLUSTY models. These latter models had Teg values spaced at roughly 500K intervals 
between 10,620K and 31,564K and a slightly larger step between 31,564K and 63,735K. The 
remaining synthetic spectra were for Kurucz models covering the Tes range 5000K-10,000K. 
In the wavelength range of interest, contributions from ^TOOOi^T or cooler T^s values are very 
small or negligible. 

Our initial study of the intermediate state adopted a variety of assumed mass transfer 
rates, together with a Standard Model accretion disk. We used the system parameters of 
Table 1, and BINSYN in a diagnostic mode, to attempt a fit to the intermediate state 
lUE spectra. With the exception of the parameters D, i, T^s^s, ^a, ^b, and H (see Table 1 
comments), the Table 1 parameters for the WD and secondary star are the same as in 
II2004. Our model divides the accretion disk into 66 annuli, with each annulus divided 
into 90 azimuthal segments. The system synthetic spectrum is the integrated sum of the 
contributions of the WD, secondary component, accretion disk, and accretion disk rim. (Note 
that, at an inclination of 7°, the accretion disk rim contributes a negligible observable effect.) 

We find that it is not possible to fit the intermediate state lUE spectra by adding a 
Standard Model accretion disk to the low state model of H2004 and using our TLUSTY 
annulus synthetic spectra. The radial temperature gradient in the Standard Model produces 
either too much fiux at short wavelengths or too little fiux at long wavelengths. In making 
this (and all subsequent) comparisons with the lUE spectra we maintained the same scaling 
divisors, for both the lUE spectra and the synthetic spectrur n, that were adopted for the 



low state comparison. Based on the measured E{B — V) = 0.0 (ISzkody &: Downes Ill982l ) we 
applied no interstellar extinction corrections to the lUE spectra extracted from the archive. 
Fig. 4 illustrates the discrepancy produced by adding an accretion disk component with 
a mass transfer rate of l.OxlO~^7Vl0yr~^. Note that this rate equals the value given by 
Hameury et al. (1988) at the top of the period gap. This mass transfer rate, and a Standard 
Model disk, produces disk T^s values ranging from 23,765K at the inner disk edge to 4222K 
at the outer disk edge, with a temperature maximum of 34,972K at r = 1.3611 xr^j. The 
absolute fiux from the accretion disk, at the Earth, is a strong function of the mass transfer 
rate and the approximate fit in Fig. 4 indicates that the adopted mass transfer rate is nearly 
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Table 1. BINSYN Model System Parameters 



parameter 


value 


parameter 


value 


M„„i(Mn,) 


0.73 


'^wd 


0.01017 


q 


0.43 


log Qy^d 


8.25 


p(d) 


0.1329 


r^(pole) 


0.28782 


^(^0) 


1.11120 


(point) 


0.40362 




98.74 


rs(side) 


0.30011 




2.74 


r^back) 


0.33258 




7.0 


log 5's(pole) 


4.95 




47,000 


log (point) 


3.61 


Teff,,(K,pole) 


2600 


log c/,(side) 


4.87 


-^wd 


1.0 


log Qsihock) 


4.68 


As 


0.5 




0.47 


bwd 


0.25 




0.01 


bs 


0.08 


H{Re) 


0.0110 



Note. — wd refers to the WD; s refers to the 
secondary star. D is the component separation of 
centers, Q is a Roche potential. Temperatures are 
polar values, A values are bolometric albedos, and 
b values are gravity-darkening exponents, r^d is in 
units of component separation, not solar radii, ra 
specifies the outer radius of the accretion disk, set 
at the tidal cut-off radius, and is the accretion 
disk inner radius, while H is the semi-height of the 
accretion disk rim. 
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correct, i.e., the mean residual over the wavelength interval plotted is small. The synthetic 
spectrum is too blue. A larger mass transfer rate increases the flux at all wavelengths, making 
the mean residual larger, and produces a synthetic spectrum that is too blue. A smaller mass 
transfer rate decreases the flux at all wavelengths, again making the mean residual larger, 
and produces a synthetic spectrum that is too blue. 

The Fig. 4 residuals suggest a departure from a Standard Model. Stationary models 
are possible (lLasotall200ll ) whose T{R) profiles differ from the Standard Model. Examples 
include systems in which irradiation by the WD is important, or syste ms in which the mass 



transfer stream directly heats the outer regions of the accre tion disk (jVisniac fc Diamond 



1989. 


1992; 


Smak 


1994 


), (but see 



flux through the outer part of the accretion disk for unspecified reasons (IKfiz &: Hubeny 



19861 ). In all of our synthetic spectra, the synthetic spectrum has absorption lines where the 
/[/i? spectrum shows emission lines. SPT95 argue that the broad emission lines in the "high" 
( V~13 — 14) state arise from the accretion disk, in contrast to narrow (Fig. 2) emission lines 
produced on the irradiated secondary when MV Lyr is in a low state. 

Two variations of the Standard Model T(R) relation affect the corres ponding system 
synth etic spectrum: (1) a change in the functional dependence of T{R) on R (jOrosz fc Wade 



20031 ) . and (2) truncatio n of the accretion disk at some inner radius larger than the WD 



radius (ILong et al. 



19941). Physically, t runcation can occur either by evaporation of the 



inner disk annuli (IHameury et al. I l2000l ) or by accretion disk interaction with a magnetic 



field associated with the WD (HL02). 

We have investigated both options. The first permits a good fit to the lUE spectra but 
the models are slightly inferior to the second option models. The second option truncates the 
accretion disk at an inner radius larger than the WD radius (ILong et al.lll994l ) and preserves 
the Standard Model T{R) functional dependence on R. Since the distance to MV Lyr is 
a fixed quantity, as the truncation radius of the accretion disk is increased, producing an 
initial reduction in flux from the model, and a resulting discrepancy with observed spectra, 
a compensating increase in the adopted mass transfer rate is necessary to restore agreement 
with the observed spectra. We started with the mass transfer rate of Fig. 4, then incremen- 
tally increased the truncation radius and made compensating increases in the mass transfer 
rate. The quality of the successive fits improved but remained perceptably discrepant (too 
blue) even up to a hmiting truncation radius of Vb/r^d = 18. At this truncation radius, 50% 
of the tidal cutoff radius, the necessary mass transfer rate becomes l.OxlO~^Al0yr~^, close 
to the one used in L1999, and one dex larger than the value given by Hameury et al. (1988) 
at the long period boundary of the period gap. We conclude that a Standard Model T(R) 
with a truncated accretion disk model produces a synthetic spectrum that is too blue and 
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1500 2000 2500 3000 
Wavelength (A) 

Fig. 4. — Fit of synthetic spectrum to intermediate state lUE spectra of MV Lyr. The 
hght hne is the lUE spectrum; the heavy hue (colored bhic in electronic edition) is the 
synthetic spectrum. The synthetic spectrum includes a Standard Model accretion disk with 
a mass transfer rate of l.Ox lO~^A^0yr~^. The synthetic spectrum was convolved with a 
Gaussian 5A FWHM broadening function to match the lUE spectral resolution. The low 
state synthetic spectrum from Fig. 1 (lower spectrum) is shown for comparison. (See the 
electronic edition of the Journal for a color version of this figure.) 
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Fig. 5. — As in Fig. 4 but with an accretion disk truncated at r = 1.7xrwd, a T(i?) corre- 
sponding to a mass transfer rate of 6.OxlO~^°A^0yr~^, and an isothermal 9500K accretion 
disk beyond r = 11.5 xr^d- The dash hne represents a corresponding black body model. 
The light line is the lUE spectrum; the heavy line (colored blue in electronic edition) is the 
synthetic spectrum. The low state synthetic spectrum from Fig. 1 (lower spectrum) is shown 
for comparison. (See the electronic edition of the Journal for a color version of this figure.) 
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so does not fit the observed MV Lyr lUE spectra. 



Our initial models maintained the WD T^s of 47,000K (H2004). ISion I ([19951) has shown 



that mass transfer can cause compressional heating of WDs in CV syst ems by 10,000K 



20,00 0K for mass transfer rates of order 1.0 x 10 ^A4qJt ^ or larger (also see lTownsley fc Bildsten 



20021 ). We tested the model sensitivity to the WD Tes by substituting a NLTE, Teg = 
52,000K, log g=8.25 WD synthetic spectrum for the 47,000K synthetic spectrum. The 
change in the system synthetic spectrum, for the wavelength range 1100A-3000A was almost 
undetectable at 3000A. In view of our much smaller mass transfer rate than those studied 
by Sion, we did not test larger Tes values and we do not consider this topic further. 

Our extensive experiments ultimately led us to a best-fitting model with a mass transfer 
rate of 6.OxlO~^°A^0yr~^, a truncation radius of = 1.7xr^^, and with T{R), which 
initially follows the Standard Model, set to 9500K for radii larger than 11.5xr„d- The T^s 
of the innermost annulus, including the correction for irradiation by the WD, is 30,195K. 
Fig. 5 shows the fit of this solution to the lUE data. The mean residual of the fit to the 
877 points between 1150A and SOOOA was 0.013 and the mean absolute residual was 0.158. 
These residuals apply to the scaled spectra. The synthetic spectrum is too blue. The dash 
line represents a corresponding black body model in which all source synthetic spectra have 
been replaced by black body spectra of the same Tes- The black body model is not blue 
enough. 



Kriz &: Hubeny I (119861 ) (their Fig. 2) provides optical thickness as a function of wave- 
length for two TLUSTY-based accretion disk models. A rough interpolation indicates that 
our models likely are optically thick. If the accretion disk remains nearly Keplerian, our 
Fig. 5 model, for equal bolometric luminosity, would correspond to a Standard Model of 
mass transfer rate higher than 6.OxlO~^°A^0yr~^, and likely near l.Ox lO~^7Vl0yr~^. 

We have not attempted to calculate irradiated models of the individual accretion disk 
annuli, a possible option in TLUSTY. Because of the disk truncation, no boundary layer has 
been included in the synthetic spectrum. 



5. Predicted change in V magnitude between the low and intermediate states 



The AAVSO light curve shown in II2004, and the Roboscope curve shown in lHoneycutt fc Kafka 



(I2OO4J ). covering 1990-2004, reveals that MV Lyr exhibits high states lasting for 2-6 months, 
followed by low states lasting for similar lengths of time. Typical high state V values 
are 12.2-12.8, while low state V values are in the range 16-18. Intermediate state mag- 
nitudes show considerable variation, and in the V range 13-15. Earlier B and V light 
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curves ( iRosino. Romano fc Marziani Ill993l . hereafter RRM) cover the interval 1968 to 1991 
and show that MV Lyr was in an extended low state from 1979 to 1989, with occasional 
transitions to intermediate states followed by returns to the low state. 

RRM tabulate B = 14.6 on Oct. 17, 1979, B = 13.9 on Nov 20, 1979, and B = 17.6 (a 
low state) on Apr 11, 1980. Since a typical B - V value for MV Lyr is -0.25 (RRM), we 
estimate V = 13.6 at the time of the lUE intermediate state observations (Dec. 2, 1979), 
roughly a magnitude below typical high state values, and at the end of the intermediate state. 
We calculate a B magnitude change, l ow state-intermediate state of AB=3.7 magnitudes. 



with AV approximately the same. The iHoneycutt &: Kafka I (120041 ) MV Lyr light curve plot 



gives a typical intermediate state-low state AV = 4.0. 

We have calculated an extension of the Fig. 5 mod el to 10,000A, and have con volved 



the spectrum with the transmission profile of the V filter (iMatthews fc Sandage Ill963l ). The 
normalized integral of this product tabulation is a measure of the theoretical V magnitude. 
This, together with a corresponding calculation with the low state spectrum, extended to 
10,000A, permits determination of the theoretical V magnitude change between the low and 
intermediate states. The Fig. 5 model produces AV = 4.91, compared with the observed 3.7 
magnitude change in B (and V) between Nov. 20, 1979 and Apr. 11, 1980 (RRM). We have 
calculated several models that fit the lUE spectra well, using elevated Tcs values in the outer 
part of the accretion disk, and all of them produce a AV that is at least one magnitude too 
large. If there is a bright spot in the intermediate state system the discrepancy is only made 
worse by adding a model of it to our calculated system model. This impasse indicates that 
a different approach is needed, and the black body plot in Fig. 5 suggests that the difficulty 
is associated with the annulus synthetic spectra. 



6. A limiting intermediate state model 

All of our models assume an outer boundary coincident with the tidal cutoff radius. Ter- 
minating the outer accretion disk boundary at a smaller value has the prospective advantage 
of reducing the flux at 5500A while having a minor effect at 3000A. We have followed this 
option to a limiting model: an isothermal 14,000K accretion disk terminated at r = 20xrwd- 
This termination radius produces a flux level that matches the observed spectrum, given 
the distance to MV Lyr. We have preserved a truncation radius at r = 1.7xr„d, although 
the flux contribution from the innermost annuli is small. The calculated AV for this model, 



low to intermediate state, is 4.06, in good agreement with the IHoneycutt fc Kafka I (120041 ) 
amplitude but larger than the RRM amplitude. Fig. 6 shows the fit of this model to the 
lUE spectrum. The required fiux level determines the accretion disk radius. The radiation 
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peak at 1300A marks the same peak of the 14,000K source synthetic spectrum of an annulus 
model. The fact that the hmiting model synthetic spectrum is too blue is a property of 
the Balmer continuum region of the source synthetic spectrum, beyond the direct control of 
BINSYN. The black body plot of Fig. 5 suggests that, if formation of the continuum takes 
place under conditions that reduce the Balmer jump and slope of the Balmer continuum, 
then the Fig. 6 model might fit the lUE spectra closely. Reducing the T^s of the isothermal 
accretion disk, to make it less blue, requires an increase in the accretion disk radius, to 
provide the necessary flux. This increases the V amplitude, producing a disagreement with 
observation. Our limiting model represents the best fit to the full array of observational 
data. 



7. Fit to the high state HST spectrum 

As part of an HST snapshot program, we obtained a 730 s exposure of MV Lyr on 
June 24, 2003 when it was in a peak high state of V~12.5. The Space Telescope Imaging 
Spectrograph (STIS) was used with grating G140L, giving a resolution near 1.2 A. Data 
reduction used the latest release of CALSTIS (V2.13b). This version takes into account the 
decaying sensitivity of the G140L grating. All the lines are in absorption (Fig. 7, discussed 
below), consistent with an optically thick disk. 

We investigated whether a Standard Model synthetic spectrum, including the option of a 
truncated accretion disk, can fit the HST spectrum. We began with the L1999 mass transfer 
rate of l.lxlO~^Al0yr~^ and no truncation and found that the predicted flux for this model 
is appreciably too large at all HST wavelengths. The same mass transfer rate and an accretion 
disk truncated at r/r^d = 6.2 leads to a synthetic spectrum that is too blue and a slightly too 
small flux level, and an accretion disk truncated at r/rwd = 5.3 produces a synthetic spectrum 
that is too blue and a shghtly too large flux level. Successively smaller mass transfer rates 
and no disk truncation produced models with flux levels still too large but with a smaller 
mean residual. At each mass transfer rate it was possible to find a pair of truncation radii for 
which the corresponding synthetic spectra bracketed the HST spectrum, one with slightly 
too large flux values, the other too small; in each case, both bracketing spectra were too 
blue. In addition to the initial test, we tested mass transfer rates of 7.0x lO^^TVlQyr^-'^, 
5.0x10-^ MqJT-\ 4.5x10-^ MqJt~\ 4:.0x10-^MqJY-\ and finally S.OxIQ-^MqJT-K In 
this latter case an accretion disk extending to the WD surface fits the HST spectrum well 
between 1150A and 1200A, but is much too blue (too little flux) longward of 1200A. Smaller 
mass transfer rates produce synthetic spectra with too small flux at all HST wavelengths. 
We conclude that Standard Model synthetic spectra, with or without truncated accretion 
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Fig. 6. — As in Fig. 5 but with an isothermal 14,000K accretion disk extending from an 
inner truncation radius at r = 1.7xr„d to r = 20.0xr„d- The hght hne is the /t/i? spectrum; 
the heavy hne (colored blue in electronic edition) is the synthetic spectrum. The low state 
synthetic spectrum from Fig. 1 (lower spectrum) is shown for comparison. (See the electronic 
edition of the Journal for a color version of this figure.) 
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disks, do not fit the HST spectrum; they are too blue. 

With the objective to flatten the theoretical spectral gradient, we selected the Standard 
Model with small errors at the short wavelength extreme, mass transfer rate of 3.0x lO^^A^oyr^^, 
and modified the T{R) relation to an isothermal outer region with elevated temperature. We 
proceeded by initially selecting an accretion disk annulus at a large radius and setting the 
Teflf values of all annuli of larger radius to the T^s of the selected annulus. The resulting 
synthetic spectrum provided an improved fit to the HST spectrum. We continued, gradually 
moving to a smaller crossover annulus and consequently higher T^g value, until we reached 
a crossover annulus of r/r^d = 14.0 and annuli Tgg = 12, OOOK at larger radii. At smaller 
radii the annulus Tcs values continued to follow a Standard Model for a mass transfer rate 
of 3.OxlO~^A^0yr^^. The spectral gradient now was about correct but the overall fiux level 
had become somewhat too large, and so we reduced the fiux level by truncating the accre- 
tion disk. Because of the high T^a value at the truncation radius, the color of the synthetic 
spectrum is mildly sensitive to the truncation radius value. We found an optimum fit with a 
truncation radius of r/r^d = 1-7 and T^s there of 43,404K. It is likely that a slightly modified 
mass transfer rate, together with slightly changed truncation radius and crossover radius and 
outer annulus Teg value would fit the HST spectrum equally well. For that reason, we do 
not attribute special significance to the fact that the truncation radius found here is the 
same as the truncation radius for the intermediate state. The large differences between the 
HST spectrum and the synthetic spectrum absorption lines lead to larger mean and mean 
absolute residuals than for the intermediate and low states. The mean residual of the fit 
to the 970 points between 1150A and 1715A was -0.44 and the mean absolute residual was 
0.71. These residuals apply to the scaled spectra. 

Figures 7-10 show the fit to the HST spectrum of a synthetic spectrum with our optimum 
high state parameters. Although the continuum is a reasonable fit, individual line fits show 
appreciable discrepancies, especially for the high excitation lines of NV, Si IV, C IV, and 
He II. Since all of our synthetic spectra assume solar composition, the discrepancies could 
arise from: (1) abundance differences from solar composition, or (2) excitation differences 
between the formation regions for the absorption lines and the continuum. 

The observed absorption lines show a slight blue shift corresponding to a velocity of 
about — 200kms~^. This effec t is apparent in the detail plots, Fig. 8 (C HI), Fig. 9 (C II, Si 



IV), and Fig. 10 (C IV, He III ISchneider. Young, fc ShectmanI fll98l[ ) and SPT95 separately 



determine a system radial velocity of MV Lyr of 7 ~ — 35kms~^ from optical wavelength 
spectra. The velocity difference (of 165km s~^) is a clear indication that the absorption lines 
are formed in a wind. Since MV Lyr is viewed nearly pole-on, this means the wind has a 
low velocity. This interpretation is consistent with the observed line broadenings as listed in 
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Wavelength (A) 

Fig. 7. — Fit of a standard model accretion disk with a mass transfer rate of 
3.OxlO"^Al0yr~^, truncated at an inner radius of l.Txrwd, and with annuh T^s values set to 
12,000K beyond r/r^d = 14.0 (thick hne, colored blue in electronic edition) to an HST high 
state spectrum of MV Lyr (hght hne). The synthetic spectrum has been convolved with a 
Gaussian 1.2A FWHM broadening function. Line identifications are in following plots. (See 
the electronic edition of the Journal for a color version of this figure.) 
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Wavelength (A) 

Fig. 8. — Detail of Fig. 7. The intermediate strength overplotted spectrum (colored red 

in electronic edition) is a synthetic spectrum of a single annulus with a Tgff of 31,564K, 
convolved with a Gaussian 1.2 A FWHM broadening function. Note the much closer fit to 
the observed Ly a (light line) than provided by the system synthetic spectrum (heavy line, 
colored blue in electronic edition). The annulus synthetic spectrum was divided by 6.5x10^ 
before overplotting, to place the synthetic spectrum over the observed spectrum. (See the 
electronic edition of the Journal for a color version of this figure.) 
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Fig. 9. — As in Fig. 8, for a different spectral region. Tfie C II and Si III labels mark 
groups of transitions. Note that the Si IV equivalent widths in the overplotted annulus 

spcctnim (31,564K; intermediate strength line, colored red in electronic edition) arc much 
greater than for the system synthetic spectrum (heavy line, colored blue in electronic edition), 
indicating that the HST spectral lines (light plotted line) are formed in a higher temperature 
environment than for the system synthetic spectrum. Note also that the observed Si IV lines 
are broadened as well as displaced to shorter wavelengths. The annulus synthetic spectrum 
was divided by 4.7x 10^ before overplotting, to place the synthetic spectrum over the observed 
spectrum. (See the electronic edition of the Journal for a color version of this figure.) 
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Fig. 10. — As in Fig. 8. Note that the equivalent widths of the C IV doublet and the He 
II line are appreciably greater in the overplotted annulus spectrum (31,564K; intermediate 
strength line, colored red in electronic edition) than in the system synthetic spectrum (heavy 
line, colored blue in electronic edition). The spectrum at the top, colored green in the 
electronic edition, represents a 38,076K annulus, shown for comparison, and with the same 
normalizing factor as the 31,564K annulus. The annulus spectra were divided by 3.9x10^ 
before overplotting, to place the 31,546K synthetic spectrum over the observed spectrum. 
When scaled to fit the observed spectrum, the 38,076K spectrum is fairly similar to the 
31,546K spectrum. Note that the C IV doublet is resolved in both annulus spectra and 
in the system synthetic spectrum. All three of these spectra have been convolved with a 
I.2A Gaussian broadening function. (See the electronic edition of the Journal for a color 
version of this figure.) 
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Table 2, e.g., the Si IV doublet near 1400A. 

Fig. 8 shows the vicinity of Ly a. Note the much better fit of the overp lotted annulus 
synthetic spectrum [T^s = 31,564K) to the observed Ly a line than is true for the system 
synthetic spectrum, whose mean Tes is appreciably lower. This indicates that the absorption 
lines are formed in a higher temperature environment (i.e., corona or wind) than is true for 
the continuum. Fig. 9 shows the Si IV doublet near 1400A ; it is too weak in the system 
synthetic spectrum, and equivalent widths for the overplotted annulus synthetic spectrum 
(Tefj = 31, 564K) are again a better fit. 

The C IV doublet (Fig. 10) in the system synthetic spectrum is also much too weak. 
The doublet is stronger in both the overplotted annulus synthetic spectra (Teff=31,564K, and 
38,076K), but not by enough when scaled to fit the observed profile. The C II line complex 
(Fig. 9) is slightly too weak in the overplotted annulus spectrum, while the equivalent width 
of the system synthetic spectrum appears comparable to the observed value. The C III line 
complex. Fig. 8 in the overplotted annulus synthetic spectrum, appears to have an equivalent 
width approximately matching the observed value, while the system synthetic spectrum is 
too weak. The varying strengths of C lines in different ionization stages and the observed 
(large) strength of N V (Fig. 8) may be due to composition effects but we cannot exclude the 
possibility th at it is an excitation effect. Large N V/C IV ratios have been seen in several 



CV systems (iGansicke et al. II2003I ) and are likely related to evolutionary effects. Note that 
H2004 determines an overall metal composition for the MV Lyr WD oi Z = 0.3x^0, but 
detailed elemental abundances, relative to solar values, of C=0.5, N=0.5, and Si=0.2. (Since 
strong excitation effects are clearly present there is no confiict between our assumed solar 
abundance for the present analysis and the WD composition determination in H2004.) 

It is of interest to compare the emission line spectrum of the intermediate state and 
the absorption line spectrum of the high state (Fig. 11). Note that emission lines in the 
intermediate state match corresponding absorption lines in the high state. As well as can be 
judged, the relative strengths of the various resonance emission lines are similar to the relative 
strengths of the same absorption lines. The HST spectrum in Fig. 11 has been displaced 
downward to superpose the C IV complex. The two broadened spectra fit accurately at the 
superposed line boundaries; if there were a radial velocity difference between the emission 
line spectrum and the absorption line spectrum as large as the displacement of the emission 
line spectrum (e.g.. Fig. 9) from the synthetic spectrum, that difference would be detectable. 
Comparable fits result when other spectral regions are superposed, although the low signal 
level of the lUE spectrum produces some ambiguity. It is of particular interest that the 
intermediate state spectrum is fiatter than the high state spectrum. The two sets of lines 
in the intermediate and high states appear to be produced in a wind with a similar velocity 



structure in the two states. 



8. Change in system luminosity, low state to high state 

The calculated V light change for our optimum model (Fig. 7), from low to high state, 
following the procedure of is 5.91 magnitudes. The observed values are; AAVSO, 6.0, and 



Honeycutt fc Kafka I (120041 ). 5.25, bracketing our calculated value. The HST observations 
provide a separate test. We used the STIS/CCD acquisition image of MV Lyr to deter- 
mine the magnitude at the time of the HST FUV observations to be 12.4 in the F28x50LP 
filter, corresponding to a monochromatic fiux (but integrated over an extended spectral re- 
gion) of 2.205±0.05xlO~^'*ergcm~^s~^A~^ at 7228.5 A. The filter has a lower transmission 
limit at 5400A, a peak response at 6OOOA, and a gradual dropoff with half-peak response 
at 8300A and an upper transmission cutoff at 10,000A. The high state model synthetic 
spectrum, when multiplied by this filter transmission and integrated, produces a calculated 
equivalent monochromatic 7228. 5A fiux, at the Earth, for the 505pc (H2004) distance to MV 
Lyr, of 2.91x10"^^ erg cm~^s^^A^^, in approximate agreement with the observed value. 



9. Discussion 



The significant result from our low state study is that, to within observational error, 
the system spectrum can be represented by a combination of a WD spectrum and an ir- 
radiated, Tefj^poie = 2600K secondary component. There is no observational evidence for 
accretion disk emission in the MV Lyr low state. Our newly calculated contribution of the 
secondary star is an interpolation among NextGen models and includes allowance for irra- 
diation of the secondary by the WD. Both L1999 and HL02 model the low state of MV 
Lyr with a truncated accretion disk that retains most of its mass from the high state, and 
emits radiation at a level consistent with residual viscosity. Our results limit possible low 
state accretion disk emission to valu e s well below those predicted by both L1999 and HL02. 
Honeycutt. Cannizzo. fc Robertson I (119941 ) discuss the occasional sudden drops in luminos- 
ity of V974 Aql and show that the drops can be explained within the accretion disk limit 
cycle theory in terms of b ack and forth propagation of cooling and he ating transition fronts. 
As with L1999 and HL02. lHoneycutt. Cannizzo. fc Robertson I (119941 ) find that the accretion 
disk undergoes little net mass loss while mass transfer from the secondary is turned off. 



HL02 discuss slow passage of an accretion disk system from an initial stable state through 
the instability region and show that outside-in outbursts would occur but are unobserved. 
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Fig. 11. — Comparison of the high state HST spectrum (upper spectrum) and the interme- 
diate state lUE spectrum (lower spectrum) of MV Lyr. The high state spectrum has been 
displaced downward by 3.0 flux units (ordinate scale) to superpose the C IV complex at 
1540A. This plot demonstrates that the strong emission hues in the lUE spectrum become 
strong absorption lines in the HST spectrum, at (as well as can be determined) the same 
wavelengths and with similar broadening. A wavelength displacement as large as that be- 
tween the HST spectrum and the synthetic spectrum in Fig. 9 would be detectable. Both 
sets of lines appear to be produced in a low velocity wind. 



-29- 



Table 2. Prominent absorption and emission lines, FWHM values 



ion 


nominal A (A) 


HST 


lUE 


low 


intermed. 


high 


C II 


1141 


10.5 




0.5 


5.8 


4.4 


C III 


1175 


6.3 


1.4 


2.6 


2.6 


3.2 


N V 


1241 


11.5 


1.2 


2.6 


1.0 


1.8 


Si II 


1260 


2.4 


1.2 


0.3 


8.0 


3.9 


Si III 


1299 


12.3 


2.0 


2.2 


1.7 


4.8 


C II 


1334 


3.0 


7.5^*^ 


1.0 


2.8 


2.7 


Si III 


1342 


9.9 




1.0 


2.7 


3.1 


Si III 


1374 


8.6 


1.5 


1.6 


1.3 


1.6 


Si IV 


1393 


6.4 


7.0**^ 


1.8 


2.0 


3.0 


Si IV 


1402 


7.4 




1.7 


2.0 


3.0 


Si III 


1501 


6.6 


7.0^ 


0.8 


2.1 


2.4 


C IV 


1549 


7.9 


10.8^ 


4.5 


3.2 


5.1 


He II 


1640 


6.6 


7.0^ 


3.4 


3.0 


3.4 



^emission line 

^this entry is the FWHM sum with the ion on the following line 



Note. — The lUE column refers to the intermediate state spec- 
trum, SWP07296. The columns designated as "low", "intermed.", 
and "high" refer to the system synthetic spectra corresponding to 
Fig. 1, Fig. 5, and Fig. 7 respectively. 
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The only solution to the VY S cl phenomenon that t hey visualize, to meet the constraint 



imposed by the case of TT Ari (iGansicke et al. Ill999l ). is a magnetic field sufficiently large 



(B>6MG) to truncate the residual low state accretion disk at a radius large enough to 
suppress accretion disk formation. However, as HL 02 note, if a magnetic field is present, 
then circular polarization should also be present. iRobinson et al.l (jl98ll ) chronicles that 
repeated attempts to measure the circular polarization in MV Lyr always found a value 
less than 0.13%. Robinson et al. also note that, if there is a magnetic field present that is 
strong enough to disrupt the accretion disk, Zeeman-splitting of absorption lines should be 
detected; it is not. Thus, there is no evidence for a strong magnetic field associated with the 
MV Lyr WD. 

Both the intermediate state /[/i? spectra and the high state HST spectra are inconsistent 
with our calculated Standard Models, having fiatter fiux profiles. Our subsequent thin disk 
simulations, using TLUSTY annulus models, differ from the Standard Model by assuming 
a ffatter T{R) profile; they are able to produce fairly good fits to the observed continua. 
Confirmation that the TLUSTY annuli represent a thin disk model is provided by the output 
data for an annulus at r/r„d = 18.5 and for a mass transfer rate of 4.0x lO~^A^0yr~^, close 
to the value of our high state model. For this annulus, whose radius is 1.45x10^'' cm, a 
Rosseland optical depth of 0.81 is located ai z = 1.94x10^ cm, one dex smaller than the 
radius. A Rosseland optical depth of 0.001 occurs aX z = 1.98x10^ cm, indicating line 
formation in a thin photosphere. 

Smak (1994) has discussed peculiar (i.e., non-Standard Mo del), TjR) distributions. The 

depa rtures from the Standard Model in the systems he cites (e.g.. lRutten. van Paradijs. fc Tinbergen 
19921 ) are in the direction of a too fiat T{R) profile, as we find for MV Lyr. Smak argues 
that the departures for systems with i > 75° result from the assumption of a fiat accretion 
disk. That explanation is not applicable for MV Lyr because of its low inclination. Smak 
suggests that, for low i cases, heating of the outer part of the accretion disk by the stream 
collision could be the explanation. Buat-Menard, Hameury & Lasota (2001) discuss the ef- 
fect of stream impact heating and tidal effects. As our high state model indicates, if stream 
impact heating is the physical cause of anomalous heating in the outer accretion disk region, 
the heating effects extend over an appreciable fraction of the accretion disk radius. This 
effect differs from a bright spot, for which the radial extent amounts to a few percent of the 
component separation (see IWarner Ill995l . sect. 2. 6. 5). We stress the importance of the low 
orbital inclination in MV Lyr and the consequential absence of eclipses. All of the systems 
cited by Si nak were sho wn to have peculiar T{R) distributions by application of the MEM 
technique (IHorne Ill993l ) and so were restricted to systems showing eclipses. The spectrum 
synthesis method used in this paper represents an alternative and independent technique 
that is applicable to both eclipsing and non-eclipsing systems. 
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The presence of high-excitation absorption hnes in the MV Lyr system imply their for- 
mation in a high temperature region above (both faces of) the accreti on disk, possibly similar 

to a chromo s phere or corona. Accretion disk coronae a re coi isidered inlMever fc Mever-Hofmeister 
Jl989l . ll994h . lLiu. Mever fc Mever-Hofmeister Nl997h . and lHameurv et al. I J200oh . We have 
presented evidence for formation of both the intermediate state emission lines and the high 
state absorption lines in a wind, apparently coextensive with the high temperature region. 
It would be of interest to searc h for a P-Cyg ni profile, possibly in the infra-red. We thank 



the re feree for this suggestion . iDrew I (119971 ) summarizes known properties of winds from 



CVs. iProga. Stone, fc Drew 



tion disks, and iPereya et al. 




discuss a model of radiation-driven winds from accre- 
discuss the stability of line-driven winds. Table 2 lists 
the FWHM wavelength ranges for the principal absorption lines in the HST spectrum. In 
most cases, several transitions contribute to a given line. The low MV Lyr i value nomi- 
nally would imply narrow disk lines, distinctly different from the observations. The observed 
broadening is consistent with lines of sight that traverse a range of radial velocities before 
reaching an optical depth of 1.0. The 200kms~^ blue shift of the absorption lines is smaller 
than typically found for winds from accretion disks. 

Our TLUSTY annuli models are in LTE and apply default values for the internal param- 
eters, leading to a negative temperature gradient perpendicular to the surface of the annulus. 
Line formation in a vertically extended wind violates our mo del assumption of hydrostatic 
equilibrium; addition of a wind model, (ILong &: Knigge II20021). rn ight provide improved rep- 
resentation of the MV Lyr spectral data. iFroning et al.l (120021 ) combined TLUSTY disk 
models with a representation of a biconical wind in a study of SS Cyg. 

We speculate that continuum formation takes place under conditions different from the 
assumptions of our annulus models and that difference may e xplain our failur e to achieve 
an accurate continuum fit, especially in the intermediate state. Hubeny ( 1989al Jbl) discusses 
conditions for hot upper layers of accretion disk annuli to exist. H e shows that the z-di r ection 
profile of viscosity strongly affects the formation of these regions. Hubeny fc Hubeny I ( 1998 ) 
discuss the treatment of viscosity as implemented in TLUSTY, in particular the steps taken 
to prevent a high temperature z-direction anomaly. The authors also emphasize the impor- 
tance of NLTE effects. It is widely be lieved that the physical c a use of acc r etion disk viscosity 
is the magnetorotational instability (IBalbus fc Hawley Ill99l[ ). iBalbus I (120021 ) discusses is- 



sues in the implementa tion of the instabilit y in accretion disk models. The magnetohydrody- 
namical simulations of IStone et al. I (119961 ) find an increase in viscosity toward the accretion 
disk surface, leading to increased z-direction heating. 

The flatter spectral profile of the intermediate state lUE spectrum, as compared to the 
HST high state spectrum (Fig. 11) is of particular interest. The isothermal accretion disk 
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extending only half way to the tidal cutoff radius in the intermediate state has the attractive 
feature that it provides a ready explanation for emission lines (the line emission region is seen 
extending beyond the edge of the accretion disk), while the same lines are seen in absorption 
in the high state (the accretion disk now extends to the tidal cutoff radius). This scenario 
indirectly supports the absence of an accretion disk in the low state, again in disagreement 
with L1999 and HL02. An explanation of the flat T(-R] profi l e, both intermediate and high , 
must differ from the case of DW UMa (jKnigge et al.l l200d : lAraujo-Betancor et al.l |2003| ) , 
also a VY Scl star. In DW UMa there is also a fiat T{R) profile, but it is explained by self- 
occultation by a puffed-up outer accretion disk rim and an orbital inclination large enough 
to produce eclipses of the WD. 



10. Summary 

(1) Based on the excellent fit of a system synthetic spectrum to FUSE, lUE, and optical 
low state spectra (H2004), any accretion disk that may be present in the low state must 
have a T^ff less than 2500K. This result agrees with a study of the CV system TT Ari 



(iGansicke et al. Ill999l ). A comparison of our results (Fig. 3) with two models (L1999 and 



HL02) of the VY Scl phenomenon demonstrates a discrepancy with those models. 

(2) Using TLUSTY annulus models, and corresponding synthetic spectra, to simulate an 
intermediate state for which lUE spectra are available, we show that a Standard Model (see 
text) cannot produce a system synthetic spectrum that satisfactorily fits the lUE spectra. A 
rough fit is possible with a Standard Model with a mass transfer rate of l.Ox lO~^A^0yr^^ 
but the model synthetic spectrum is too blue. A fairly good fit is possible with a mass 
transfer rate of 6.x lO~^°7W0yr~^, a truncation radius of r = l.Txr^d, and with r(i?)=9500K 
beyond r = ILSxr^d- But this model produces a change in V magnitude between the low 
and intermediate states that is about one magnitude larger than the observational value. 
An isothermal 14,000K accretion disk with a radius that is half of the tidal cutoff radius 
produces agreement with the V magnitude change; this model is slightly bluer than the lUE 
spectra. 

(3) Study of a high state demonstrates that no combination of mass transfer rate and 
truncation radius can produce a Standard Model system synthetic spectrum whose contin- 
uum satisfactorily fits our recent HST spectrum. All synthetic spectra are too blue. A 
fairly good fit is possible for a mass transfer rate of 3.0x lO^^A^oyr"^, a truncation radius of 
r = 1.7xr„d, and T(i?)=12,000K beyond r = 14.0xrwd- This model produces a calculated 
system flux at 7228. SA that is in reasonable agreement with a measurement of HST optical 
flux. 
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(4) The absorption lines in our MV Lyr HST spectrum show a range of excitations, and 
corresponding temperatures, that are higher than the nearly flat continuum would imply. 
This result is consistent with line formation in a large z high temperature region. Anomalous 
line strengths of C, N, Si, and He have a possible interpretation in terms of composition 
effects, but excitation effects cannot be disentangled. 

(5) The slight blue shift of the absorption lines (~ — 200kms~^), as compared with the 
system radial velocity of — 35kms~^, indicates that the absorption lines are formed in a low 
velocity wind, which is coextensive with the large z high temperature region. There is no 
apparent difference in velocity between the emission line intermediate state spectrum and 
the high state absorption line spectrum; the wind and high temperature structures appear 
to be present in both the intermediate and high states. 

(6) The difference between the physical conditions under which the absorption and emis- 
sion lines form in a wind and the physical conditions assumed by our hydrostatic equilibrium 
annulus models, used to calculate the synthetic spectra, may explain our difficulty in fitting 
the high state and intermediate state continua accurately. The same problem may explain 
the failure of our synthetic spectra to represent the high excitation absorption lines in the 
high state HST spectra. 
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